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1. Introduction 

The theoretical calculation of the hadronic B — > nlv form factor f+(q 2 ) is a key ingredient in 
the determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \V u b\ from B — > nlv 
exclusive semi-leptonic decay. The value of \V u b\ characterizes the strength of the quark-flavor 
changing b — >■ u transition, and can be obtained by combining the hadronic B — > nlv form factor 
f+(q 2 ) with experimental measurements of the differential decay rate via 

dT(B^nlv) G 2 F \V ub \ 2 f 2 2 2 2 2 2 i3/2, , 2 , |2 n 
= 192^4 [(, " s + ^"° " 4Wfl ^ (L1) 

where the momentum transfer q^ = pg — and we neglect the mass of the outgoing lepton. The 
form factor f+(q 2 ) encodes nonperturbative QCD dynamics, and can only be computed precisely 
from first principles using lattice QCD. Indeed there have been two 2+1 flavor lattice calculations 
of f + (q 2 ) done by the HPQCD [1] and FNAL/MILC Collaborations [2]. Both groups use the MILC 
gauge configurations. 

The precise calculation of \V u b\ constrains the apex of the CKM unitarity triangle. There is 
a persistent puzzle, however, between independent determinations of |V M &|. First, the exclusive 
determination from B — > nlv and the inclusive determination from B — > XJv (where X u is any 
charmless hadronic final state) differ at the level of more than 3 a [3]. Second, most experimental 
measurements of BR (B — > TV) combined with lattice-QCD input for /g give a higher value of 
\V u b\ than both |V M fc|exci and IV^Iinci [4, 5], although the recent Belle measurement of BR (B — > TV) 
is lower and more in line with the Standard-Model prediction, albeit still with large errors [6]. 
Given this situation, an independent lattice calculation using a a different gauge action is desired to 
address this puzzle. 

2. Methodology 

The B — > nlv hadronic weak matrix element is parameterized by the form factors f+{q 2 ) and 

Mq 2 ) as 

W\B) = f + (q 2 ) U +£- ?*^V) /o(<7 2 )^V, (2.1) 

where the b — > u vector current "V^ = iuy^b. We calculate numerically the form factors f\\ and f±, 
which are more convenient for the lattice calculation: 

(n\^\B) = ^2rrr B [v^f\\{E n )+p^f L {E n )], (2.2) 

where = Pg/ms and p±= p% — {Pn • In the B- meson rest frame, these form factors are 
proportional to the hadronic matrix elements of the temporal and spatial vector current: 

(n\Y°\B) {n\r\B) 1 

ill = — /~ — > 71 = — ik — —r- y 1 ^) 

The desired form factor /+ can be obtained by following relation: 

f + (q 2 ) = -j^==[f\\(E n ) + (rriB — E„)f±(E„)]. (2.4) 
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Table 1: Lattice simulation parameters. The results reported here use the ensembles specified in boldface 
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Figure 1: Three-point correlation functions for computing the B — > 7r/v form factor (left) and renormaliza- 
tion factor Z v b (right). The single and double lines correspond to light- and b-quark propagators, respectively. 
The spectator light quark is labeled I 1 and the daughter light quark is labeled /. Shaded cricles denote the 
gauge invariant Gaussian smeared source/sink for the b-quarks. 

We match the lattice amplitude to the continuum matrix element using the mostly nonpertur- 
bative method of Ref . [7] : 

W\B)=Zl{K\V»\B), Z»=p»Jzfz». (2.5) 

The flavor-conserving renormalization factors Z v h and Zy are computed nonperturbatively on the 
lattice and the factor p is computed at one loop in mean-field improved lattice perturbation the- 
ory [8]. Most of the heavy-light current renormalization factor comes from Z v b and Zy, such that 
p is expected to be close to unity [9]. 

We improve the b — > u vector current through &(a,s a )- At this order we need only compute 
one additional matrix element with a single-derivative operator. We calculate the improvement 
coefficient at 1-loop in mean-field improved lattice perturbation theory. We have not yet included 
the G (a) -improvement term in the results shown in these proceedings. 

3. Computational setup 

Our computation of the B — > %l v form factor is performed on 2 + 1 flavor domain- wall fermion 
and Iwasaki gauge-field ensembles generated by the RBC and UKQCD Collaborations with sev- 
eral values of the light dynamical quark mass at two lattice spacings, the coarser a ps 0.11 fm 
{a- 1 ps 1.73 GeV) and the finer a ps 0.086 fm (a" 1 ps 2.28 GeV) [10, 11]. The set of gauge-field 
configurations used for this project is summarized in Table 1. The results reported here use the 
ensembles specified in boldface. On the finer ensembles we compute two quark propagators on 
each configuration with their temporal source locations separated by J/2 to increase the statistics. 
Periodic boundary conditions are imposed in the time direction. In this proceedings, we present 
first results obtained from the coarser ensembles and a subset of the finer ensembles. 
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For the light quarks, domain-wall valence quark propagators are generated on each config- 
urations with several partially quenched masses in order to enable good control over the chiral 
extrapolation. For the bottom quark, we use the relativistic heavy quark (RHQ) action [12] to 
remove the large discretization error introduced by the large bottom quark mass [13]. The RHQ 
action is given by 



Srhq = £ 9n \ mo + 7oA) - — ^ + % 



y-D 



a{Df 



(3.1) 



where tuning the bare-quark mass m^a, the clover coefficient cp, and the anisotropy parameter % for 
the b quark is needed [12, 14]. Here we employ values determined nonperturbatively in Ref. [15]. 

Figure 1 shows the the three-point correlation functions needed in this project. The left dia- 
gram of Fig. 1 depicts the computation of the B — > nlv form factor. The source of the pion (with 
both zero and nonzero momenta) is located at the origin. The B-meson is at rest at t sl ±. We compute 
this three-point function with a unitary spectator mass and multiple partially quenched daughter- 
quark masses that enable interpolation to the physical strange-quark mass and extrapolation to the 
physical average up-down quark mass. The right diagram of Fig. 1 is used to obtain the renormal- 
ization factor Z v h . We fix the the spectator mass to be m q = m s in order to reduce the statistical 
eiTors because Z v b is independent of the light spectator-quark mass. For Zy, we use the value ob- 
tained by the RBC/UKQCD Collaborations in Ref. [11]. To reduce excited-state contamination, 
we employ a gauge-invariant Gaussian smeared (sequential) source for the B-meson in all lattice 
correlators. In order to reduce the statistical error of the three point function, we calculate the three 
point function with the B-meson at t^k and (T — t SI ±) and average the results. 

4. Two-point and three-point fits 

The pion energy and B-meson mass are obtained from the following two point functions: 

C$(t,p) = £^-(^(?,x)^(0,0)), (4.1) 

.V 

cf (0 = £(AM^(o,o)) 5 (4.2) 



where (J n and O'b are interpolating operators for the pion and B-meson. We compute the effective 
energies as 



E e ff{p) = cosh 



-l 



C 2 {t,p)+C 2 (t + 2,p) 



(4.3) 



C 2 (t + l,p) 

The left plot of Fig. 2 shows the effective pion energies for momenta through (pL/2n) 2 = 3 on the 
the coarser ami = 0.005 ensemble. Similarly, the right plot of Fig. 2 shows the effective B-meson 
mass calculated with a Gaussian smeared source and point sink. 

In the lattice-QCD simulations, the form factors /jj at and /J 31 can be obtained from the follow- 
ing ratios of correlation functions: 



R 



B— >7l 



(Msnk) 



Cf (t)C B (t snk -t) 



2En 



e -E n t e -m B t 



(4.4) 
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Figure 2: Effective energy/mass plots for the pion (left) and B-meson (right) on the coarser ami — 0.005 
ensemble. Shaded bands show the fit results with jackknife statistical errors and fit ranges. For the pion, data 
points for four spatial momenta (pL/2n) 2 = 0, 1,2,3 are indicated by different colors/symbols. 
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Figure 3: The ratio R^ K (left) and R^ K (right) with f snk = 20 on the coarser ami 
ranges and fit results with jackknife statistical errors are shown as horizontal bands. 



0.005 ensemble. Fit 



4 at = lim R%^(t,t SJ± ), (4.5) 
/?= lim ±-R?^(t,t SBk ), (4.6) 

0«f«f s nk P' K 



where 



is the vector current matrix element on the lattice. 
Fig. 3 shows the results for R^ K and R^ 71 calculated with a source-sink separation t snk = 20 
on the coarser ami = 0.005 ensemble. For R^ K , we average over equivalent spatial momenta. 
Plots of these ratios for the other ensembles and partially quenched masses look similar. We deter- 
mine _/jj at and /J 31 by fitting to a plateau in the region <C t <C t snk where we expect the excited-state 
contributions to be negligible and obtain a good correlated # 2 /d.o.f. and p-value. 
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HZj? = 10.037(34), p = 12% 
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Figure 4: Determination of the renormalization factor Z y b . with t sl ± — 20 on the coarser am/ = 0.005 
ensemble. The shaded band shows the fit result with jackknife statistical error. 

5. Renormalization factor Z bb 

Values of Zy b can be computed nonperturbatively using the charge-normalization condition 
Zffi (B s \V bb '°\B s ) = 1 where V bb '° is the b -> b lattice vector current. We obtain Z v b by fitting the 
following ratio of correlators to a plateau in a region <C t <C t SI ±: 

Zy b (t,t sn k) = C^'^snk) I Cf^ Bs (t,tsr±). (5.1) 

The determination of Z bb on the coarser ami = 0.005 ensemble is shown in Fig. 5. At tree level in 
mean-field improved lattice perturbation theory, Z bh is given by [12, 13] 

Z v b = Mo exp(M 1 ), M\ = log [1 +mo] , m Q = — - (1 + 3£)(1 - — ). (5.2) 

Using mo = 8.45, £ = 3.1 and mo = 0.8757 for the coarser ami = 0.005 ensemble [11, 15], we 
obtain a tree level estimate Z v b = 10.606, which is in good agreement with the nonperturbative 
determination. 

6. Form-factor results and outlook 

After multiplying the results for R^ 71 and R^f n with the renormalization factor Z bl , we ob- 
tain the form factors fu and f± shown in Fig. 5. (For now the p-factor is set to unity.) After adding 
the p-factors and 6(a) -improvement, we will extrapolate to the physical quark masses and con- 
tinuum and interpolate in E\ using chiral perturbation theory. Calculations on the remaining finer 
ensembles are underway. 

We will then extrapolate our synthetic form-factor data over the full kinematic range (down 
to q 2 = 0) using the model-independent z-parameterization [16-19]. Our result will provide an 
important independent check on existing calculations using staggered light quarks. 
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Figure 5: Form factors fn and fx in coarse lattice units 024. The red, bue and green symbols indicate data 
on the coarser am/ = 0.005, ami = 0.01 and finer ami = 0.004 ensembles, respectively. Different symbols 
indicate different partially quenched masses. 
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